Multifractality in stock indexes: Fact or fiction? 
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Abstract 

Multifractal analysis and extensive statistical tests are performed upon intraday minutely 
data within individual trading days for four stock market indexes (including HSI, SZSC, 
S&P500, and NASDAQ) to check whether the indexes (instead of the returns) possess 
multifractality. We find that the mass exponent r(q) is linear and the singularity a(q) is 
close to 1 for all trading days and all indexes. Furthermore, we find strong evidence showing 
that the scaling behaviors of the original data sets cannot be distinguished from those of 
the shuffled time series. Hence, the so-called multifractality in the intraday stock market 
indexes is merely an illusion. 
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1 Introduction 



Econophysics is an emerging interdisciplinary field applying concepts, theories, 
and tools borrowed from statistical physics, nonlinear sciences, applied mathemat- 
ics, and complexity sciences to understand the complex self-organizing behaviors 
of financial markets id AS Hi- This field has become to flourish since the pioneer- 
ing work of Mantegna and Stanley on the scaling behavior in the dynamics of the 
Standard & Poor's 500 index [Q], which is closely related to the Pareto-Levy law 
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proposed by Mandelbrot in the description of cotton price fluctuations flg]. Econo- 
physicists have uncovered remarkable similarities between financial markets and 
turbulent flows yj,|4j]. Such analogues include (but not limited to) the evolution of 
probability densities of financial returns 10] based on the variational theory in turbu- 
lence flBESEH], inverse statistics in stock markets I 12 L 13 L 14]motivated by the 
inverse structure function analysis of velocity II 151 llq Il7l Il8l 1911 . scale-invariant 
distribution of multipliers defined from volatility of equities Il20ll and from dissi- 
pating energy [|2lL l22l. l23L 12411. and intermittency and multifractality of asset returns 

HQ. 



Indeed, the multifractal nature of equity returns is one of the most important styl- 
ized facts. A small part of this literature contains the studies on the foreign exchange 
rate flESHHHSH^ldprice H], commodityprice Q, returns 



of stock price or indexes 11321 . 133L l34l 135 



38L [39L |40|, |4J], and so on. We 



note that the quantity price (or its logarithm) in financial markets is the analogue 
of velocity in turbulence. Similarly, the counterpart of velocity difference in fluid 
mechanics is the asset return. In this framework, it is natural that numerous multi- 
fractal analyses have been carried out on the returns for financial equities similar to 
the velocity differences for turbulent flows. 



However, there are exceptions, where analysis is performed on several indexes di- 
rectly rather than their variations (the returns) and the presence of multifractality in 
the several indexes is claimed Il42ll43l . 14411 . Specifically, they performed multifractal 
analysis on the intraday high-frequency data of Hang Seng Index (HSI), Shanghai 
Stock Exchange Composite Index (SSEC), and Shenzhen Stock Exchange Com- 
posite Index (SZEC) within individual trading days. The extracted "multifractal" 
spectra f(a) were then utilized to predict abnormal price movements and serve as 
a risk measure in risk management. It seems to us that a careful scrutiny on the 
obtained multifractality should be undertaken based on the extremely narrow spec- 
tra of the singularity a. Two problems arise, casting doubts on the aforementioned 
analysis ll4~5ll . 



Firstly, based on the multifractal theory, there exists a constant a(t) for each mo- 
ment t such that the investigated measure fx on the neighbor B(t, I) of x scale with 
I when the scale / — > 0, 

fi(B(t,i))~rv. (i) 

The measure /i is singular at arbitrary moment t with the singularity strength being 
a(t). When /x is defined as the sum of index prices within a given time interval, 
/x (B(t, I)) is approximately proportional to /, that is, a(t) m 1 for all t. This sug- 
gests that the measure fx does not possess multifractal nature. This inference is fur- 
ther supported by the fact that the span of singularity strength Aq = ct max — a 
in the real data 0,00]. 



mm ~ 



Secondly, in the analysis of multifractality in turbulence or high-frequency financial 
data, the moment order q should not be greater than 8 in order to make the partition 
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function converge. Specifically, it is shown that the size of a time series should be 
no less than one million to ensure the estimate of its eighth order partition function 



statistically significant 11191 146H. The situation is similar for high-frequency financial 
data [20]. Hence, it is of little significance to compute partition function for higher 
orders. In the analysis of minutely (or five-minute) data within a time period of one 



day 11421 143L 14411. the size of the intraday high-frequency data is no more than 240 
while the moment order is taken to be —120 ^ q ^ 120. This usually broad interval 
of q casts further doubts on the reported multifractality in the indexes. 

Despite of the specific considerations discussed above, it is worthwhile to put fur- 
ther comments in general on the investigation of multifractality in financial data. 
The multifractal features in financial series have attracted great interests, however, 
the origin and significance of the extracted "multifractality" is less concerned. On 
one hand, it has been shown that an exact monofractal financial model can lead to 



an artificial multifractal behavior [47]. On the other hand, a time series of the price 



fluctuations possessing multifractal nature usually has either fat tails in the distri 



bution or long-range temporal correlation or both B48I1 . However, possessing long 



memory is not sufficient for the presence of multifractality and one has to have a 



nonlinear process with long-memory in order to have multifractality I149ll . In many 



cases, the null hypothesis that the reported multifractal nature is stemmed from the 



large fluctuations of prices cannot be rejected 1150m 



In this work, we focus on the presence of multifractal feature in stock market in- 
dexes and testing whether the obtained empirical multifractality stems from random 
fluctuations. To address these issues, we adopt the bootstrap approach by shuffling 
the intraday index series and perform multifractal analysis on them. The results are 
compared with that from original data. This paper is organized as follows. In Sec.[2l 
we describe the data sets we investigate. The basic multifractal method is explained 
in detail in Sec.|3] Multifractal analysis of the data sets is presented in Sec. 51 Sta- 
tistical bootstrapping tests are conducted in Sec.[5l Finally, Sec.|6]concludes. 



2 Data sets 



To gain a more profound insight into the multifractality in intraday stock market 
indexes, we investigate four important indexes, i.e., the Hang Seng Index (HSI), 
Shenzhen Stock Exchange Composite Index (SZSC), Standard & Poor's 500 Index 
(S&P 500), and the National Association of Securities Dealers Automated Quota- 
tion (NASDAQ). HSI and SZSC are selected since they were used in the original 



work of this topic 11421. 143L 14411 . Both the Hongkong Stock Exchange and Shen- 
zhen Stock Exchange are emerging markets. The S&P 500 and NASDAQ that are 
representative of mature stock markets are chosen for comparison. 

The data have been recorded at each minute in trading days. The HSI index covers 



3 



from Jan. 2, 1997 to May 28, 1997, the SZSC index is from Nov. 12, 2001 to Aug. 
17, 2006, the S&P 500 index is recorded from Jan. 2, 1997 to Feb. 26, 1999, and 
the NASDAQ index ranges from Aug. 18, 2000 to Oct. 30, 2000. Eliminating the 
weekend, holidays and, the days having recording errors, there are 101 days for the 
HSI data, 1 149 days for the SZSC data, 448 days for the S&P 500 data, and 45 days 
for the NASDAQ data, respectively. 



3 Method 



We use the box counting method following the work of H42I . I43L 14411 to investigate 
the multifractal nature of the index series of each trading day. Denote the intraday 
index series as {I(t) : t = 1, 2, • ••, T}, where T = 240 for HSI and SZSC, T = 405 
for S&P 500, and T = 390 for NASDAQ, respectively. For a given box size /, we 
obtain N = T/l boxes and construct a measure [i on each box as follows, 

i 

fx(n; I) = fj, {[{n - 1)1 + 1, nl\) = £ J[(n - 1)/ + z] , (2) 

i=l 

where [(n-l)l+l, nl] is the n-th box and I E [1, 2, 3, 4, 6, 10, 15, 20, 30, 40, 60, 80, 
120, 240] for HSI and SZSC, / e [1, 3, 5, 9, 15, 27, 45, 81, 135, 405] for S&P 500, 
and / E [1, 2, 3, 5, 10, 13, 15, 26, 30, 39, 78, 130, 195, 390] for NASDAQ, respec- 
tively. The sizes / of the boxes are chosen such that the number of boxes of each 
size is an integer to cover the whole time series. 

We then construct the partition function x q as 

/*(n;i) 



x q (i) = E 



n=l 



E m =iMra;0 



(3) 



and expect it to scale as 

x Q (i) ~ r (q) , (4) 

which defines the exponent r(q). The local singularity exponent a of the measure 



jj, and its spectrum f(a) are related to r(q) through a Legendre transformation 115 ill 

| a = dT(q)/dq 
\f(a) = qa-r(q) 



In order to keep the comparability of our results with those in |44j], we also pose 
-120 < q ^ 120. 

When fj,(n; /i(m; I) < 1 and q ^> 1, the estimate of the partition function 
X will be very difficult since the value is so small that it is out of the memory. 
To overcome this problem, we can calculate the logarithm of the partition function 
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lnx g (i) rather than the partition function itself. A simple manipulation results in 
the following formula, 



N 



n=l 



li{n;l) 



max{/i(m; I)} 



gin 



max{/i(m; /)} 

m 



(6) 



where max{/x(m; Z)} is the maximum of /i(m; /) for m = 1, 2, • • • , iV. 



4 Multifractal analysis 



Four dates (Jan. 8, 1997 for HSI, Nov. 26, 2001 for SZSC, Feb. 10, 1997 for S&P 
500, and Aug. 22, 2000 for NASDAQ) are taken as examples to show the results of 
multifractal analysis. Figure Q] shows the dependence of the partition function x?(0 
on the box size I for different values of q in log-log coordinates. Excellent power- 
law scaling of Xg(0 with respect to / has been observed and the scaling range covers 
all the selected values of /. The solid lines are the best linear fits to the data. 
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Fig. 1. Plots of XqQ) as a function of the box size I for different values of q in log-log 
coordinates. The solid lines are the least-squares fits to the data using linear regression (in 
log-log coordinates) corresponding to power laws, (a) HSI, (b) SZSC, (c) S&P 500, and (d) 
NASDAQ. 

The scaling exponents r(q) are given by the slopes of the linear fits to lnx q (l) 
with respect to In/ for different values of q. Figure [2] plots the dependence of 
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the mass exponents r(q) as a function of the moment order q. One observes that 
there is an evident linear relationship between r(q) and q for all the four exam- 
ples. The solid lines are the least-squares fits to the data. The slopes of the lines 
are respectively a = 1.000 ± 0.001 for HSI, a = 1.000000 ± 0.000003 for SZSC, 
a = 1.00000 ± 0.00001 for S&P 500, and a = 1.0001 ± 0.0001 for NASDAQ, re- 
spectively. All the corresponding correlation coefficients of the linear fits are equal 
to 1.0000. Furthermore, the linear relationships are also hold for other trading days. 
Therefore, there is no evidence of nonlinearity in the functions r and the intraday 
stock market index do not exhibit multifractal nature. Since a(q) = dr(q)/dq, we 
expect that a(q) « 1 for all q, as expected in our discussion in Sec. CD 




Fig. 2. Dependence of the scaling exponent r{q) on the order q. The solid lines are the 
least-squares fits to the data, (a) HSI, (b) SZSC, (c) S&P 500, and (d) NASDAQ. 



Figure |3] presents the multifractal singularity spectra f(a) obtained through Leg- 
endre transformation of r(q) defined by Eq. ©. The curves in Fig. |3] have the 
geometrical features of the conservable multifractal spectra [flU^l], which makes 
them look as if there is sound evidence for the presence of multifractality. However, 
when looking at the disperseness of the sigularity strength Aa — a max — a m i n , we 
find that Aa is very close to zero. It is well-known that Aa is an important pa- 
rameter qualifying the width of the extracted multifractal spectrum. The larger is 
the Aa, the stronger is the multifractality. According to Fig. [3l even in the case 
of -120 ^ g ^ 120, Aa < 0.002 for NASDAQ. One can see that the values of 
Aa for other indexes are much smaller than that of NASDAQ. This observation 
indicates that there is no multifractality in stock market indexes. 
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Fig. 3. (Color online) Multifractal spectra f(a) obtained by the Legendre transform of r(q) 
for different indexes. 

5 Statistical tests for multifractality 
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Fig. 4. (color online) Comparison of multifractal spectra extracted from real and shuffled 
stock marker indexes. The solid lines are the real data, while the dotted lines are the shuffled 
data, (a) HSI, (b) SZSC, (c) S&P 500, and (d) NASDAQ. 

We access further the statistical significance of the empirical multifractality in the 
sprit of bootstrapping. For a given intraday time series, we reshuffle the series to 
remove any potential temporal correlation and carry out the same multifractal anal- 
ysis as for the original data. For the four examples discussed in Sec. H we compute 
the multifractal spectra of ten reshuffled time series for each index. The results are 
illustrated in Fig. HI where the solid lines are associated with real stock market 
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indexes, while the dotted lines are obtained from the shuffled data of the corre- 
sponding indexes. We find that the multifractal spectra of the real indexes f(a) and 
that of the shuffled data / rn d(a r nd) are almost overlapping together in Figj4](b) and 
(c). Although the solid lines and the dotted line can be distinguished clearly in Figj4] 
(a) and (d), the differences between a and a rn d are ignorable. In other words, the 
multifractal nature in the real indexes is insignificant in these examples. 

For each intraday time series, we shuffle the data for 1000 times. The associated 
multifractal spectra are obtained. For each singularity spectrum, we calculate two 
characteristic quantity, Aa and F = [f(a min ) + /(a max )]/2. Figure [5] shows the 
scatter plots of F rn( j for the shuffled data versus the corresponding Aa rn d for the 
four example trading days. Clear linear relationship between F rn d and Aa rn d for 
each case is observed and we have 



rnd 



A;Aa rnd + b 



(7) 



where k = -30.31 and b = 1.00 for HSI, k = -30.10 and b = 1.00 for SZSC, 
k = -30.05 and b = 1.05 for S&P 500, and k = -30.31 and b = 1.06 for 
NASDAQ, respectively. The open circle in each plot of Fig. [5] presents the values 
of F and Aa for the real data. 
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Fig. 5. Scatter plots of the dependence of the shuffled F rn d and the corresponding Aa m d. 
(a) HSI, (b) SZSC, (c) S&P 500, and (d) NASDAQ. 

Two striking facts emerge from Fig. [5] First, the 4000 points of (Aa rnd , -F rn d) col- 
lapse on a same linear line since the values of k and b are identical for the four 
plots. Second, the four points of (Aa, F) for the four real data sets also locate on 
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the same line. For other trading days, we have observed similar phenomena, which 
put further evidence on our conclusion that the real and reshuffled time series have 
undistinguishable scaling behaviors. 

The values of Aa and F for each original time series are compared with the aver- 
ages (Aa rn d) and (F rn d) of the 1000 corresponding shuffled data sets. The results 
for the four indexes are illustrated in Fig. [6l The solid line is the main diagonal 
y = x. We find that Aa ~ (Aaw) and F ~ (-F m d) for all cases, which implied 
that the multifractal spectra of the shuffled data are very close to that of the real data 
and the f(a) curves of real index data can be completely interpreted by the random 
fluctuations of the original data sets. We stress that there are no extreme values in 
the intraday index prices so that one can not attribute the observed multifractality 
to tail fatness that is absent in the present case. Hence, the multifractal property in 
high-frequency stock market indexes obtained by partition function method is not 
statistically significant. It is just an illusion. 
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Fig. 6. Comparison of F and Aa obtained from the shuffled data and the real data, (a) HSI, 
(b) SZSC, (c) S&P 500, and (d) NASDAQ. 

In the presence case, to test the presence of multifractality amounts to testing 
whether the local singularity exponent a ^ 1, or Aa ^ 0. As a last step, we 
impose a very strict null hypothesis to investigate whether the f(a) spectrum is 
wider than those produced by chance. The null hypothesis is the following: 

H :Aa^ Aa rnd . (8) 

We can compute the p-value, which is the probability that the null hypothesis is 
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true. The smaller the p- value, the stronger the evidence against the null hypothe- 
sis and favors the alternative hypothesis that the presence of of multifractality is 
statistically significant. The false probability is estimated by 

p x = Pr(Aa ^ A« rnd ) . (9) 

Under the conventional significance level of 0.05, the multifractal phenomenon is 
statistically significant if and only if pi ^ 0.05. While pi > 0.05, the null hypothe- 
sis cannot be rejected. A similar null hypothesis can be described as follows: 

H : F > F Tnd , (10) 

where the false probability is 

p 2 = Pr(F ^ F rnd ) . (11) 

Using the conventional significance level of 0.05, the multifractal phenomenon is 
statistically significant if and only if p 2 ^ 0.05. 

For the four examples shown in Fig. [5l we find that pi — 1 and p 2 — 1 for HSI, 
p x = 0.108 and p 2 = 0.109 for SZSC, p x = 0.452 and p 2 = 0.456 for S&P 500, and 
Pi = and p 2 = for NASDAQ. Obviously, we can not distinguish the real data 
from the shuffled data beside NASDAQ for the chosen trading days. We also find 
that pi ps p 2 for all the trading days. More generally, Table Q] shows the statistical 
tests for the all the each trading days. About half of the trading days can not pass 
the statistical inference, indicating that multifractality is absent in the those trading 
series. 
Table 1 

Statistical tests for the presence of multifractal nature in the four indexes investigated. 



Indexes 


HSI 


SZSC 


S&P 500 


NASDAQ 


Percentage of p\ < 0.05 


54.6% 


56.1% 


54.4% 


53.9% 


Percentage of p 2 ^ 0.05 


54.4% 


55.8% 


53.6% 


53.6% 



6 Conclusion 

We have investigated the multifractal features in intraday minutely high-frequency 
stock market indexes (including HSI, SZSC, S&P 500, and NASDAQ) for indi- 
vidual trading days. The resultant scaling functions r(q) have been confirmed to 
be linear and the singularities a are close to 1 so that Aa is close to 0. This 
analysis implies that there is no multifractality in the indexes. Further evidence 
based on bootstrapping technique shows that that the scaling behavior of the shuf- 
fled data is undistinguishable from that of the raw data. Specifically, we find that, 
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(1) almost all points (Aq, F) of the raw data sets locate on the same straight line 
Fmd — — 30Aa rnd + 1 extracted from the points (Aa rnd , F rnd ) of the shuffled data; 

(2) for each time series, Aa « (Aa rnd ) and F ~ {F rni \)\ and (3) the two rather 
strict null hypotheses cannot be rejected for about half of the time series. There is 
thus no doubt that the reported multifractal nature in the indexes of HSI and SZSC 



11421 143L 14411 is not a fact but a fiction. This conclusion is further verified by two 
indexes (S&P 500 and NASDAQ) in a developed stock market. We believe that our 
analysis and conclusion apply for other market indexes or common stock prices 
when one concerns intraday stock prices or indexes rather than their returns. 

In addition, we cast doubts on the efforts to use this illusionary multifractal feature 
to forecast the stock market [430 and to define a risk index for risk management 



114411 . However, to be more conservative, we do not deny the potential usefulness 
of those techniques proposed based on some nonexistent properties. The idea to 
use multifractal nature to predict or to manage risks in stock markets should be 
investigated based on the returns or other alternative financial quantities. After all, 
one cannot build a palace on a sand beach. 
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